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The density-of-states (DOS) and one-particle spectral function A(k, w) of the one- and two-orbital 
models for manganites, the latter with Jahn- Teller phonons, are evaluated using Monte Carlo 
techniques. Unexpectedly robust pseudogap (PG) features were found at low- and intermediate- 
temperatures, particularly at or near regimes where phase-separation occurs as T— >0. The PG 
follows the chemical potential and it is caused by the formation of ferromagnetic metallic clusters in 
an insulating background. It is argued that PG formation should be generic of mixed-phase regimes. 
The results are in good agreement with recent photoemission experiments for Lai,2Sri.8Mn207. 

PACS numbers: 71.10.-w, 75.10.-b, 75.30.Kz 
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Manganese oxides have been intensively studied in re- 
cent years due to their colossal magnetoresistance (CMR) 
effect 0, where small changes in magnetic field induce 
large changes in their resistivity pdc- The CMR behavior 
may originate in the nontrivial interplay of the charge, 
spin, orbital, and lattice degrees of freedom in these ma- 
terials. Among the theories proposed for the CMR ef- 
fect are those based on the recently observed intrinsic 
tendencies of models for manganites toward mixed-phase 
formation In this context, clusters of one phase em- 
bedded into the other, typically involving ferromagnetic 
(FM) metallic and FM or antiferromagnetic (AF) insu- 
lating phases, are believed to substantially increase pdc 
and the compressibility k of the system Q . The effect oc- 
curs even in stable metallic regimes at low-temperatures, 
as long as the density is close to those where phase- 
separation occurs Q. A variety of experimental results 
are in qualitative agreement with these ideas and evi- 
dence is accumulating that manganites have important 
microscopic inhomogeneities, particularly in regimes of 
relevance for the CMR effect gj| . 

Recently, high-energy resolution angle-resolved pho- 
toemission (ARPES) measurements for the bilayer com- 
pound Lai.2Sri.8Mn207 (hole density x=0.4) have been 
reported The most remarkable result was the exis- 
tence of a pseudogap at the chemical potential /i, sub- 
stantially larger than the analog feature reported for the 
cuprates. The main purpose of this paper is to analyze 
whether similar PG characteristics appear in the DOS of 
manganite models when analyzed using unbiased compu- 
tational many-body methods [^j . The results are surpris- 
ingly robust, showing that indeed a prominent pseudo- 
gap exists in the DOS, as long as the regime explored has 
mixed-phase characteristics. The conclusions are general 
and they should apply to other compounds with similar 
microscopic phase-separation tendencies as well. 

The one-orbital Hamiltonian is given by H = Hde + 



Haf , with the Double-Exchange (DE) term defined as 
Hde 
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where cl^ creates an Cg-fermion at site i with spin cr, and 
Si is the total t2g-spin, assumed localized and classical 
(|Si|=l). Jh>0 is the Hund coupling, and the rest of the 
notation is standard. Haf = J'X](ij) ' ^ direct 
AF coupling between localized spins. The two-orbitals 
Hamiltonian is Hkjt = Hk + Hjt + Haf • The first term 
containing the kinetic energy and Hund coupling is 



H, 
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where a,b=l,2 are the two eg-orbitals, u are unit vec- 
tors along the main axes, u labels those axes, Si = 
X)aa/3 '^ia.a'^af3Ciiif3 IS the Spin of the mobilc fermions, and 
the rest of the notation is standard In this paper, the 
unit of energy will be t=l in Eq.(l), tf^ = 1 in Eq.(2a), 
and J' wiU be fixed to 0.05 §. The couphng with JT- 
phonons (QP=-Qp=Qp', and QP^QP^qI')) is 



Hjt = A X cLQr'cib. 
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Classical phonons are assumed for simplicity. The eg- 
fermionic density (n) is adjusted with /i, for both models 
Eqs.(l,2) (the hole-density is x = 1 — (n)). To study 
their properties, a Monte Carlo (MC) algorithm for the 
classical spins and phonons and an exact diagonalization 
of a one-electron problem for a given spin-phonon back- 
ground were used. The method allows dynamical calcu- 
lations in real-time without uncontrolled analytical con- 
tinuations, and it has been described in detail in previous 
publications The DOS N(a;) was obtained adding the 
one-particle spectral functions A(k, w), and it is directly 
related with the DOS of photoemission experiments. 

Fig. la contains ^{uj) of the one-dimensional (ID) one- 
orbital model Eq.(l) at large Hund coupling (Jh=oo for 
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simplicity) and (n)~0.88 (x ^ 0.12). In this regime, MC 
calculations of (n) vs n at zero temperature (T) similar 
to those reported in Refs. |||,^ revealed the presence of 
phase-separation (k^oo) for 0.77<(n)<1.0. As a conse- 
quence, the finite-T state investigated in Fig. la presents 
a large compressibility. Size effects are small, and the re- 
sults are expected to be representative of the bulk limit. 
The most noticeable feature of Fig. la is the deep min- 
imum at Lo—fi, which clearly develops as T is reduced. 
Similar behavior has been observed for other parameters 
(Jh, J', T and (n)) investigated here, as long as k is large. 
Fig. lb shows two-dimensional (2D) results, still for the 
one-orbital model. Here phase-separation at T=0 was 
found in the interval 0.86<(n}<1.0. Fixing T to a low 
value, and varying (n), PG behavior is once again ob- 
served where T=0 phase-separation occurs. In Figs.la-b 
the PG is always centered at /i, i.e. the DOS is not rigid 
but adjusts with (n) and T. Also note that at densities 
with marginally stable T=0 ground-states, PG remnants 
are observed (insets of Figs.la-b), as long as k remains 
large. This occurs in a narrow (n)-range near phase- 
separated regimes. 
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FIG. 1. Pseudogap formation in the one-orbital model, as 
a function of (n) and T. (a) Density of states N(lj) vs cu-fj, at 
Jh=oo and ID. The overall density is (n)~0.88 (x ~ 0.12). 
Starting from the top at u}-fi=0, the temperatures T of the 
first four lines are 1/10, 1/17, 1/25, 1/40, and the lattice 
size is L=30. The dotted line with crosses was obtained at 
T=l/40, but with L=1QQ showing the absence of strong size 
effects. The inset contains the DOS for L=100, T=l/40, and 
{n)~0.77, a marginally stable density as T— >0; (b) N(cj) vs 
uj-jj, in 2D at Jh=cx3, using a 10 x 10 cluster and T=l/30. The 
four lines from the top at u-fi^Q correspond to (n) ~0.90, 
0.92, 0.94, and 0.97. The inset has resuhs for {n)=0.86, a 
marginally stable density at T=0. 

Pseudogap features appear also in the two-orbital 
model with JT-phonons, again as long as k is robust. For 
instance, working at large Jh and A, Fig. 2a shows results 
at (n)~0.7, varying T. For the parameters of Fig. 2a, MC 
studies showed that T=0 phase-separation occurs in the 
interval 0.5<(n)<1.0 between spin FM-phases that differ 



in the orbital-order pattern (uniform vs staggered) 
Fig. 2a illustrates the PG development as T is reduced, 
similarly as for the one-orbital model. An analogous phe- 
nomenon also occurs at low-density (Fig. 2b) where T=0 
phase-separation between spin FM- and AF-phases was 
numerically observed for 0.0<(n)<0.3 1^. At low (n), re- 
sults without electron-phonon coupling (Fig. 3a) are simi- 
lar to those with strong lattice coupling A = 1.5. Vestiges 
of the PG also appear at (n) ~ 0.4 (Fig. 2c), which has 
a stable T=0 ground-state based on previous studies [|| . 
A similar situation occurs at A = 1 and (n) 0.7, near 
the phase-separation boundary (Fig. 3b). 
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FIG. 2. Pseudogap behavior of the ID two-orbitals model 
at large electron-phonon coupling, as a function of (n) and T. 
(a) N(tj) vs uj-iJL with Jh=8, A = 1.5, L=20, and (n) ~ 0.7 
(x ~ 0.3). The thin-dashed, thick-dashed, and solid lines cor- 
respond to T=l/5, 1/10 and 1/20, respectively; (b) Same as 
(a) but at Jh = oo, L=22 and (n) ~ 0.2. The thin-dashed, 
thick-dashed, and solid lines correspond to T=l/5, 1/10 and 
1/25, respectively; (c) Same as (a) but for T=l/15 and 
(n) ~ 0.4, a stable density as T ^ 0. 



To investigate the origin of the PG, in Fig. 3c the 
charge-charge correlation in momentum-space N(q) is 
shown in the region of interest. At low hole-density, N(q) 
is enhanced at small-q in the PG regime, indicating the 
existence of extended charged structures. Similar conclu- 
sions are reached by noticing that the average FM-cluster 
size increases as T is reduced (Fig. 3d). The overall study 
of the data accumulated here lead us to conclude that the 
PG formation is correlated with large charge fluctuations, 
which occur both with and without phonons in regions 
that phase-separate as T— >0, and in their vicinity. 

Although the relation between PG formation and 
mixed-phase tendencies is clear, the sharpness of the PG 
feature and its location following as (n) is varied re- 
quires additional analysis. For this purpose, several t2g- 
spin MC-configurations (snapshots) of the ID one-orbital 
model (Fig. la) were individually studied (the reasoning 
and results described below apply with minor changes 
to the two-orbital model with JT-phonons and in any 
dimension). The analysis of the one-electron energies re- 
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vealed that most of the MC spin configurations produce 
the PG feature, and, thus, a single MC-generated snap- 
shot should be enough to understand the phenomenon. 
Visual analysis indicated that the configurations with PG 
contain FM clusters of an intermediate size between 2 to 
6 lattice spacings, immersed in a background with AF- 
correlations (short-ranged since here T is comparable to 
J'). To understand why this spin arrangement produces 
a PG, consider for simplicity the low-(n) limit. Since at 
Jh = oo, the high and low density limits are identical, 
while at finite but large Jh they are very similar, our 
analysis applies to (n) close to 1 simply changing elec- 
trons by holes. A cartoon-like representation of a typical 
MC t2g-spin snapshot is in Fig. 4a 0. The electrons are 
correlated with the classical spins such that they mostly 
populate the FM regions (Fig. 4b), and, thus, aligned t2g- 
spins provide to the electrons an attractive effective po- 
tential (Fig. 4c). When extra electrons are added, FM 
clusters (with one or more electrons) are created, and 
new (occupied) levels appear below /x. Then, the one- 
particle spectrum at w < for one MC classical spin 
configuration contains energies corresponding to quasi- 
localized electrons in each FM-well, with a finite band- 
width ( "cluster" band) caused by the different shapes of 
those wells, plus tunneling among them. The AF-regions 
contribute to w > /x. As a consequence, the resulting 
DOS (Fig.4d) contains a pseudogap at uj = fj,. 
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FIG. 3. (a) N{uj) vs ui-jj, for the ID two-orbitals model 

without phonons (A = 0) at low-density (n) ~ 0.1, Jh = oo, 
and T=l/20. Shown are sizes L=40 and 22, denoted by solid 
and dashed lines, respectively; (b) Same as (a) but at A = 1.0, 
T=l/25 and (n) ~ 0.7. This is at the boundary of the T=0 
phase-separated region [3]. The solid and dashed lines cor- 
respond to L=40 and 50, respectively; (c) N(q) for the ID 
one-orbital model with Jh=oo, T=1/30, and L=30. Open 
circles, full diamonds, open squares, full triangles, and open 
triangles, denote (n) = 0.75, 0.84, 0.87, 0.90, and 0.94, respec- 
tively; (d) Size of typical FM clusters for the ID one-orbital 
model. Sp(j) (normalized to Sp(0)) is the classical spin-spin 
correlation (Si ■ Si+j) (j > i) measured only when Sf(1) is 
positive as a way to isolate FM clusters. 



The present study also included the individual one- 
particle spectral functions A(k, w). Figs.5a-b show the 
2D one-orbital model at large Jh, low-T, and low hole- 
density. The data displays the PG feature, which ap- 
pears at all values of the momenta which have appre- 
ciable weight at uj—fj,. Similar behavior is obtained for 
the two-orbitals model at large A (Fig. 5c). These re- 
sults are compatible with the explanation of PG behav- 
ior based on mixed-phase tendencies, since the clusters 
formed in such a state are not distributed forming a reg- 
ular pattern. Then, there is no preferred momentum for 
the PG to occur, unlike in a charge-density-wave state. 
The results of Fig. 5 are in excellent qualitative agree- 
ment with ARPES measurements for Lai.2Sri.8Mn207 
where the PG was observed all along the "Fermi sur- 
face". Note also the dispersive character of the results 
of Fig. 5, in qualitative agreement with the ARPES data 
They originate mainly from the AF-regions that in- 
duce an effective DE-hopping toff =t(cos(0/2)) which for 
the temperatures in Figs.5a-b is approximately t/2. 
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FIG. 4. Schematic explanation of the PG formation in the 
low-(n) limit and for the one-orbital model. In (a) a typical 
MC t2g-spin configuration is sketched. In (b) , the correspond- 
ing Cg-density is shown, with electrons mostly located in the 
FM regions. In (c), the effective potential felt by the electrons 
is presented. A (populated) "cluster band" is formed (thick 
line). In (d) the resulting DOS is shown. 

Other unexpected effects in ARPES data are also re- 
produced in Fig. 5. For instance, the width of the peaks 
found in Ref. Q were as large as 0.5-1.0 eV, similar to 
the widths in Fig. 5 which are of order l-2t, with t esti- 
mated to be 0.2-0.4 eV |^]. The large width is caused by 
the disordering influence of localized spins and phonons, 
particularly within the AF-clusters which have only par- 
tially ordered spins at the T's analyzed here. The peaks 
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do not sharpen as is reached, as in experiments , and 
the concept of quasiparticle seems not apphcable in the 
models studied here. 

Motivated by the similarities of our results with those 
of ARPES studies for the bilayer manganites here it 
is conjectured that the x=0.4 state of Lai.2Sri.8Mn207 
should present microscopic inhomogeneities, both above 
and below Tc. This is compatible with a recently 
reported neutron-scattering-based phase diagram [ p^ , 
where coexistence of FM (metal) and AF (insulator) fea- 
tures were observed at x=0.4 and low-T. Our predic- 
tion is also compatible with other neutron scattering re- 
sults that reported mixed AF/FM characteristics and 
short-range charge ordering due to cluster formation, for 
the same compound above Tc. 
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FIG. 5. A(k,(j;) at Jh=oo. (a) is for a 2D one-orbital 
case with T=l/30, (n) ~ 0.92 (x ~ 0.08), and on a 12 x 12 
cluster with periodic boundary conditions (PBC). Shown are 
results along (0, 0) to (tt, 0); (b) Same as (a) but along (0, 0) to 
(tt, tt); (c) Results for the 2D two-orbitals model with T=l/10, 
A=1.5, (n) ~ 0.7 (x ~ 0.3), on a 10 x 10 cluster, along (0,0) 
to (7r,0), with PBC. 



The general character of our explanation of the com- 
putational results suggests that PG behavior should be 
present in other compounds if they are in mixed-phase 
regimes. Previous studies [|[|^ suggest that this may oc- 
cur in the small x region of Lai-xSrxMnOa at low- and 
intermediate-T, as well as within the FM metallic phase 
at low-T close to the metal-insulator transition where n 
is expected to be large Q. On the other hand, at x^O.3 
this compound behaves more as a regular metal and it 
should not have a PG. In general, where Tc is the largest 
is where PG tendencies should be the weakest, correlated 
with weak CMR effects. In Lai_xCaxMn03 at low-T, the 
interesting region is again small x, but also x ~ 0.5 where 
phase-separation has been observed . If the insulating 
character of this compound above Tc is caused by micro- 
scopic inhomogeneities, PG features should also appear 
in a wide range of compositions at intermediate-T. 

Summarizing, simple models for manganites studied 
computationally have a DOS with robust pseudogap fea- 



tures in regions of parameter space with large compress- 
ibility, particularly low hole-doping. N(/x) is drastically 
reduced as T is lowered in this regime. The PG was found 
for all momenta along the "Fermi surface" , in agreement 
with experiments. The effect is caused by the forma- 
tion of FM metallic clusters in an insulating background, 
precursors of the T=0 phase-separation tendencies which 
occurs both with and without JT-phonons. Mixed-phase 
tendencies are crucial for PG formation and its associated 
Pdc increase. This conclusion seems independent of the 
details of the models, and may occur even if cluster for- 
mation is caused by 1 /r-Coulomb interactions on phase- 
separated regions. The close relation between N(^) and 
Pdc expected here differs drastically from theories based 
on Anderson localization, or in a simple hopping reduc- 
tion by spin fluctuations t^t(cos(0/2)), where N(/x) is 
approximately unchanged . The mechanism discussed 
in this paper can also be operative in non Jahn- Teller 
compounds with tendencies to phase-separation, such as 
diluted magnetic semiconductors. 
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